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Abstract 

SARAH is a Mathematica package for studying supersymmetric models. It calculates for a given model the 
masses, tadpole equations and all vertices at tree-level. This information can be used by SARAH to write 
model files for CalcHep/CompHep or FeynArts/FormCalc. In addition, the second version of SARAH can 
derive the renormalization group equations for the gauge couplings, parameters of the superpotential and 
soft-breaking parameters at one- and two-loop level. Furthermore, it calculates the one-loop self-energies 
and the one-loop corrections to the tadpoles. SARAH can handle all N = 1 SUSY models whose gauge sector 
is a direct product of SU(N) and 17(1) gauge groups. The particle content of the model can be an arbitrary 
number of chiral superfields transforming as any irreducible representation with respect to the gauge groups. 
To implement a new model, the user has just to define the gauge sector, the particle, the superpotential and 
the field rotations to mass eigenstates. 
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1. Introduction 



The LHC has started its work and will hopefully find soon the first clues for physics beyond the stan- 
dard model (SM). Supersymmetry (SUSY) is the most prominent and well studied extension of the SM. 
However, most studies were done in the context of the Minimal Supersymmetric Standard Model (MSSM) 
[lHfjj]. Therefore, many computer tools can handle the MSSM out of the box but it demands some effort to 
implement new models. Also the analytical expressions for all possible interactions, renormalization group 
equations or self-energies for the MSSM are nicely presented in literature (see e.g. Of course, this 

is not the case for all possible extensions of the MSSM. The Mathematica package SARAH is supposed to 
close this gap. The first version already allowed a comprehensive analysis of supersymmetric models by 
calculating the tree-level masses and tadpole equations as well as all interactions of the model. The obtained 
information can either be written in I^TgX files or used to generate model files for CalcHep/CompHep 
and FeynArts/FormCalc [13[. To calculate all of these results, only the minimal amount of information 
about a model is needed: the gauge sector, the particle content, the superpotential and the field rotations. 
The new version of SARAH goes one step further. First, the set of possible models which can be handled has 
been significantly extended. SARAH is no longer restricted to chiral superfields in the fundamental represen- 
tation, but can work with any irreducible representation of SU(N). Second, SARAH provides now functions 
for the calculation of the one-loop masses and the renormalization group equations (RGEs) at one- and two- 
loop level: it calculates the anomalous dimensions for all chiral superfields and the /3-functions for all gauge 
couplings, superpotential parameters, soft-breaking parameters and vacuum expectation values (VEVs). 
Furthermore, it calculates the one-loop self-energies of all fields as well as for the one-loop corrections to 
the tadpoles. This information can be easily used to get the radiative corrections to the masses at one-loop 
level. 

Before we discuss the new features of SARAH 2, we give a brief, general introduction to SARAH: in sec. [5] we 
explain the installation and the first evaluation of a model. Afterwards, in sec. [3] we show the main features 
for obtaining tree-level results and producing model files for CalcHep/CompHep or FeynArts/FormCalc. In 
sec. HI we discuss the possibilities to derive one- and two-loop RGEs, before we show in sec. [5] how to cal- 
culate the one-loop corrections to one- and two-point functions. Finally, we explain how to implement new 
models in SARAH in sec. [6] 

The appendix contains supplementary information about the check for gauge anomalies ( |Appendix A[ ), the 
calculation of the Lagrangian in component fields ( |Appendix B[ ) and the derivation of gauge group factors 
( Appendix C[ ). In addition, we list all formulas used for the calculation of the RGEs and self-energies in 
|Appendix D show our conventions for the MSSM in |Appendix E and |Appcndix F| before we summarize 
the changes in comparison to the first version of SARAH in Appendix G| 



2. Download, installation and first evaluation 

SARAH can be downloaded from 

\protect\vrule widthOpt\protect\href {http : //theorie . physik . uni-wuerzburg . de/\string~f nstaub/sarafi . htm! 

The package should be extracted to the application directory of Mathematica. This directory under Linux 
in 

home/user/ . Mathematica/ Applications/ 

and 

Mathematica-Directory\AddOns\Applications\ 

in Windows. 

Initially, the package itself consists of two directories: the directory Package contains all Mathematica 
package files, while in the directory Models the definitions of the different models are located. During the 
work, a third directory called Output is generated by SARAH. The results of different calculations as well as 
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the model files for the diagram calculators are written to this directory. A fourth directory LaTeX contains 
IATj^X packages which are needed for the appropriate output. 

A comprehensive manual (sarah.pdf) is included in the package archive and can also be found on the web 



page and on arXiv [12] . In addition, a file (models.pdf) with information about all models delivered with 
the package is part of the archive. Furthermore, a file with a short introduction to the main commands is 
included (Readme.txt) as well as an example for the usage (Example .nb). 
After the installation, the package is loaded in Mathematica via 

« " sarah-2 . 0/SARAH . m" 

and a supersymmetric model is initialized by 
Start ["Modelname"] ; 

Here, Modelname is the name of the corresponding model file, e.g. for the minimal supersymmetric standard 
model the command would read 

Start ["MSSM"] ; 

or for the next-to-minimal supersymmetric model in CKM basis 
Start ["NMSSM","CKM"] ; 

is used. In the following, we refer for all given examples the model file of the MSSM. Our conventions 
concerning the fields definitions and rotations in the MSSM are given in |Appcndix E| Even if the meaning 
of most symbols used in the examples should be intuitively clear, we list the internal names for parameters 
and particles in |Appendix F In addition, we give in |Appcndix G an overview about the main changes 



happened in the second version of SARAH in comparison to the first version presented in |11| 



3. Tree-level calculations 

When a model is initialized using the Start command, this model is first checked for gauge anomalies 
and charge conservation. If not all checks are fulfilled, a warning is printed. More information about the 
different checks is given in |Appcndix A| Afterwards, the calculation of the complete Lagrangian at tree- 
level starts. The performed steps are presented in Appendix B When this calculation is finished, several 



tree-level results can easily be obtained. 

3.1. Particle content 

To get an overview of all particles of the different eigenstates, 

Particles [Eigenstates] 

is used. This can be, for instance, Particles [GaugeES] or Particles [EWSB] for the gauge eigenstates or 
the eigenstates after electroweak symmetry breaking (EWSB), respectively. The output is a list with the 
following information about each particle: (i) name of the particle, (ii) type of the particle (F for fermion, S 
for scalar, V for vector boson, G for ghosts, A for auxiliary field), (iii) number of first generation, (iv) number 
of last generation, (v) indices of the particle. Fermions are listed using Weyl and not Dirac spinors. For 
instance, the entry for the gauge eigenstates of the left-down quark reads 

{FdL, 1, 3, F, {generation, color}} 
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3.2. Masses and tadpole equations 

The masses and tadpole equations are derived automatically during the evaluation of a model. In this 
regard, the masses or the entries of a mass matrix are calculated as second derivative of the Lagrangian 

with respect to the considered fields & . The tadpoles 7$ are the first derivative of the scalar potential with 
respect to the different VEVs Vi 

dV 

The user has access to both information by using the command MassMatrix [Particle] for the mass matrix 
of Particle and TadpoleEquation[VEV] for the tadpole equation corresponding to the vacuum expectation 
value VEV. 

Examples. 

1. Higgs mass matrix. The (1, 2)-entry of the mass matrix of the scalar Higgs in the MSSM is saved 
in MassMatrix [hh] [[1,2]]. This returns 

-(gl~2*vd*vu)/4 - (g2~2*vd*vu)/4 - B[\[Mu]]/2 -conj [B [\ [Mu] ] ] /2 

2. Squark mass matrix. In the same way, the (l,l)-cntry of the 6x6 down squark mass matrix is 
obtained by MassMatrix [Sd] [[1,1]]. The output is 

(-3*g2~2*(vd~2 + vu~2) + gl~2*(vu~2 - vd~2) + 24*mq2[l,l] + 
12*vd~2*sum[jl, 1, 3, conj [Yd[jl, l]]*Yd[jl, l]])/24 

3. Squark mass matrix with generation indices as variable. To get the result for the 2x2 down 
squark matrix without the explicit insertion of generation indices, MassMatrixUnexpanded [Sd] [[1,1]] 
is used. The output is 

(Delta [cml , cnl] * (- ( (gl~2+3*g2~2) * (vd~2-vu~2) *Delta [gml ,gnl] ) 

+ vd~2*sum[jl,l,3,conj [Yd [j 1 ,gnl] ] *Yd [j 1 ,gml] ] + 12* (2*mq2 [gml ,gnl] )))/24 

4. Tadpole equation. The tadpole equation corresponding to = is obtained by TadpoleEquation [vd] 

and reads 

(8*mHd2*vd + gl~2*vd~3 + g2~2*vd~3 - gl~2*vd*vu~2 - g2~2*vd*vu~2 - 4*vu*B [\ [Mu] ] + 
( 8*vd*\ [Mu] *conj [\ [Mu] ] - 4*vu*conj [B [\ [Mu] ] ] ) /8 == 



3.3. Vertices 

The vertices are calculated as partial derivatives of the Lagrangian with respect to the involved fields 
and applying afterwards the vacuum conditions. The vertices can be calculated in two ways. Either it is 
possible to calculate the vertices for a specific choice of external particles or to calculate all vertices of the 
complete model at once. The former task is evolved by 

Vertex [{Particles} , Options] ; 

The argument of this function is a list with the external particles. The options define the set of eigenstates 
(Eigenstates — > name) and usage of relations among the parameters (UseDependences — > True/False). 
In the result, the coefficients corresponding to different Lorentz structures are separately listed. If possible, 
the expressions are simplified by using the unitarity of rotation matrices, the properties of generators and, 
if defined, simplifying assumptions about parameters. 
All vertices for a set of eigenstates are calculated at once by 



MakeVertexList [Eigenstates , Options] ; 
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This searches for all possible interactions present in the Lagrangian and creates lists for the generic sub- 
classes of interactions, e.g. VertexList [FFS] or VertexList [SSVV] for all two-fermion-one-scalar in- 
teractions and all two-scalar-two-vector-boson interactions, respectively. If effective theories are consid- 
ered, six-particle interaction (SixParticlelnteractions — > False) or all higher-dimensional operators 
(ef f ectiveOperators — > False) can be switched off during this calculation. Those interactions slow down 
the computation and they are sometimes not needed. 

Examples. 

1. One possible Lorentz structure. Vertex [{hh, Ah, Z}] leads to the vertex of scalar and a pseudo 
scalar Higgs with a Z-boson 

{{hh[{gtl}], Ah[{gt2>], VZ[{lt3}]}, 
{ ( (ZA [gt2 , 1] *ZH [gt 1 ,1] - ZA [gt2 , 2] *ZH [gt 1 , 2] ) * (g2*Cos [ThetaW] +gl*Sin [ThetaW] ) ) /2 , 
Mom[Ah[{gt2}] , lt3] - Mom[hh[{gtl}] ,lt3] » 

The output is divided in two parts. First, the involved particles are given, second, the value of the 
vertex is given. This second part is again split in two parts: one is the Lorentz independent part and 
the second part defines the transformation under the Lorentz group. 

2. Several possible Lorentz structures. Vertex [{bar [Fd] ,Fd,hh}] is the interaction between two 
d-quarks and a Higgs: 

{{bar[Fd[{gtl, ctl}]], Fd[{gt2, ct2>] , hh[{gt3>]}, 
{ ( ( -I ) *Delt a [ct 1 , ct2] *Delt a [gt 1 , gt2] *ZH [gt3 , 2] *Yd [gt2 , gt 1] ) /Sqrt [2] , PL} , 
{ ( ( -I ) *Delt a [ct 1 , ct2] *Delt a [gt 1 , gt2] *ZH [gt3 , 2] *Yd [gt 1 , gt2] ) /Sqrt [2] , PR}} 

Obviously, there are three parts: one for the involved particles and two for the different Lorentz 
structures. PL and PR are the polarization projectors Pj, = |(1 — 75), Pr = i(l + 75). 

3. Changing the considered eigenstates and using Weyl fermions It is also possible to calculate 
the vertices for Weyl fermions and/or to consider the gauge eigenstates. For instance, 

Vertex [{fB, FdL, conj[SdL]}, Eigenstates -> GaugeES] 

gives 

{{fB, FdL[{gt2, ct2}] , conj [SdL[{gt3, ct3}]]}, 
{((-I/3)*gl*Delta[ct2, ct3] *Delta[gt2, gt3] ) /Sqrt [2] , 1}} 

4. Using dependences. With Vertex [{conj [Se] , Se, VP}, UseDependences -> True] g\ and gi 

are replaced by the electric charge e. This and similar relations can be defined in the parameter file 
(see sec. EH]) . 

{{conj [Se[{gtl}]] , Se[{gt2}], VP[{lt3}]}, 

{(-I)*e*Delta[gtl,gt2] , -Mom [conj [Se[{gtl}]] , lt3] +Mom [Se [{gt2}] ,lt3]}} 

5. Fixing the generations. It is possible to give the indices of the particles already as input 
Vertex [{hh [{1}] , hh[{l}], Ah [{2}], Ah [{2}]}] 

leads to 

{{hh[{l}], hh[{l}], Ah [{2}], Ah [{2}]}, 
{(-I/4)*(gl~2 + g2~2) *Cos[2*\ [Alpha] ]*Cos[2*\ [Beta] ] , 1}} 

Obviously, the given definition of the mixing matrices for the Higgs fields was automatically inserted. 
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3.4- Output for diagram calculators and RTfiK 

CalcHep/CompHep. CalcHep and CompHep are well known and widely used programs for calculating cross 
sections and decay widths via a diagrammatic approach at tree-level. The model files produced by SARAH 
support both Feynman gauge and unitarity gauge. Furthermore, SARAH can split interactions between four 
colored particles in a way that they can be handled by CalcHep/CompHep and also models with CP violation 
are possible. The model files for CalcHep/CompHep are created by 

MakeCHep [Options] ; 

The options define, whether the Feynman gauge should be included (FeynmanGauge — > True/False) and 
if CP violation should be considered (CPViolation — > True). Also the splitting of specific four-scalar 
interactions can be suppressed as long as they are not colored (NoSplitting — > list of fields). In addition, 
the running of the strong coupling constant can be included as it is usually done in the standard CalcHep 
files (UseRunningCoupling — > True). 

FeynArts/FormCalc. FeynArts is a Mathematica package for creating Feynman diagrams and the corre- 
sponding amplitudes. This information is afterwards used by FormCalc to simplify the amplitudes and 
square them by using FORM. In contrast to CalcHep/CompHep FeynArts/FormCalc can deal also with loop 
diagrams. Beside the standard model file for FeynArts, SARAH writes a second file including supplementary 
information about the model: all defined dependences, the numerical values for the parameters and masses, 
if they are available, and special abbreviations to speed up the calculations with FormCalc. The model files 
are generated via 

MakeFeynArts ; 

WTft&. The generated P/L3X files include all information about a model for one set of eigenstatcs: particle 
content, mixing matrices, tadpole equations, RGEs, one-loop self-energies and one-loop corrections to the 
tadpoles as well as all interactions. The I^TjtX output using the standard functions of Mathematica is not 
really readable when dealing with long expressions. Therefore, special functions were developed for a better 
typesetting. For each vertex, the corresponding Feynman diagram is automatically drawn using FeynMF flEl ]. 
In addition, a batch file for Linux and for Windows are written by SARAH to simplify the compilation of the 
different F>TeX files including all Feynman diagrams. The command for producing the P/TgX output is 

MakeTeX [Options] ; 

One option is to disable the output of the Feynman diagrams (FeynmanDiagrams — > False), another to use 
a shorter notation for the interactions (ShortForm — > True). 

4. Renormalization Group Equations 

SARAH calculates the RGEs for the parameters of the superpotential, the soft-breaking terms, the gauge 
couplings at one- and two-loop level and the VEVs. This is done by using the generic formulas of (il. Il7l - fig | 
which we have summarized in | Appendix D.l| 

The calculation of the RGEs can be started after the initialization of a model via 
CalcRGEs [Options] ; 

Options. The options offer a possibility to disable the calculation of the two-loop RGEs (TwoLoop — > 
False). Another option is to handle the number of generations of specific chiral superfields as variable 
(VariableGenerations — > list of fields). This might be used for models which include chiral superfields 
much heavier than the SUSY scale to make the dependence on these fields obvious. Normally, the /3-function 
are written in a compact form using matrix multiplication, as explained below. This can be switched off by 
the option NoMatrixMultiplication — > True. The last option offers the possibility to read the results of 
former calculations (ReadLists — > True) 
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GUT normalization. The gauge couplings of abelian gauge groups are often normalized at the GUT scale. 
Therefore, it is possible to define for each U(l) gauge coupling the GUT-normalization by the corresponding 
entry in the parameters file. See sec. 16.2.21 for more information. 

Results. The RGEs are saved in different arrays in Mathematica whose names are given in brackets: anoma- 
lous dimensions of all superfields (Gij), trilinear superpotential parameters (BetaYijk), bilinear superpoten- 
tial parameters (BetaMuij), linear superpotential parameters (BetaLi), trilinear soft-breaking parameters 
(BetaTijk), bilinear soft-breaking parameters (BetaBij), linear soft-breaking parameters (BetaLSi), scalar 
soft-breaking masses (Betam2ij), gaugino soft-breaking masses (BetaMi), gauge couplings (BetaGauge) and 
VEVs (BetaVEVs). 

All entries of these arrays are three-dimensional: the first entry gives the name of the parameter, the second 
one the one-loop /3-function and the third one the two-loop /3-function. Furthermore, the results for the 
RGEs of the scalar masses are simplified by using abbreviations for often appearing traces (see also Q). 
The definition of the traces is saved in the array TraceAbbr. In TraceAbbr [ [1] ] all abbreviations needed 
for the one-loop RGEs are given, and in TraceAbbr [ [2] ] those are for the two-loop. 
The results are also saved as text files in the directory 

PackageDirectory/Output/Modelname/RGEs/ 

Matrix Multiplication. Generally, the results contain sums over the generation indices of the particles in 
the loop. SARAH always tries to write them as matrix multiplications, in order to shorten the expressions. 
Therefore, new symbols are introduced: 

• MatMul [A,B,C, . . .] [i, j] : (ABC...)ij. Matrix multiplication, also used for vector-matrix and 
vector-vector multiplication. 

• trace [A,B,C, . . .] : Tr(ABC . . .). Trace of a matrix or product of matrices. 

• Adj [M] : Aft. Adjoint of a matrix. 

• Tp[M] : M T . Transposed of a matrix. 

As already mentioned, the usage of matrix multiplication can be switched off with the corresponding option. 
In addition, it is automatically switched off, if the model contains a parameter with three generation indices. 

Examples. 

1. /3-function of Yukawa couplings. The Yukawa couplings of the MSSM are saved in BetaYijk. The 

first entry consists of 

BetaYijk [[1,1]] : Ye[il,i2] , 

i.e. this entry contains the /3-functions for the electron Yukawa coupling. The corresponding one-loop 
/3-function is 

BetaYijk [[1,2]] : 

(-9*gl~2*Ye [il , i2] )/5-3*g2~2*Ye [il , i2] +3*trace [Yd, Adj [Yd] ] *Ye [il , i2] + 
trace[Ye,Adj [Ye]]*Ye[il, i2] +3*MatMul [Ye , Adj [Ye] , Ye] [il , i2] 

The two-loop /3-function is saved in BetaYijk [ [1 , 3] ] but we skip it here because of its length. 

2. /3-function of soft-breaking masses and abbreviations for traces. The first entry of Betam2ij 
corresponds to the soft-breaking mass of the selectron 

Betam2ij [[1,1]] : me2[il,i2] 

and the one-loop /3-function is saved in Betam2i j [[1,2]]: 
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(-24*gl~2*MassB*conj [MassB] +10*gl~2*Trl [1] ) *Kronecker [il , i2] /5 + 
4*mHd2*MatMul[Ye,Adj [Ye]] [il , i2] +4*MatMul [T [Ye] ,Adj [T [Ye] ] ] [il,i2] + 
2*MatMul[me2,Ye,Adj [Ye]] [il , i2] +4*MatMul [Ye , ml2, Adj [Ye] ] [il , i2] + 
2*MatMul[Ye,Adj [Ye] ,me2] [il,i2] 

The definition of the element Trl [1] is saved in TraceAbbr [[1,1]]: 

{Trl[l], -mHd2 + mHu2 + trace [md2] + trace [me2] - trace [ml2] + 
trace [mq2] - 2*trace [mu2] } 

3. Number of generations as variable. With 
CalcRGEs [VariableGenerations -> {q}] 

the number of generations of the left-quark superfield is handled as variable. Therefore, the one-loop 
/3-function of the hypercharge couplings reads 

(63*gl~3)/10 + (gl~3*NumberGenerations [q] )/10 

4. No matrix multiplication. Using matrix multiplication can be switched off by 
CalcRGEs [NoMatrixMultiplication -> True] 

The one-loop /3-function for the electron Yukawa coupling is now written as 

sum [ j 2 , 1 , 3 , sum [ j 1 , 1 , 3 , conj [Yd [ j 2 , j 1] ] *Yu [i 1 , j 1] ] *Yd [ j 2 , i2] ] + 
2*sum [ j 2 , 1 , 3 , sum [ j 1 , 1 , 3 , con j [Yu [ j 1 , j 2] ] *Yu [ j 1 , i2] ] *Yu [i 1 , j 2] ] + 

sum [ j 2 , 1 , 3 , sum [ j 1 , 1 , 3 , con j [Yu [ j 2 , j 1] ] *Yu [i 1 , j 1] ] *Yu [ j 2 , i2] ] + 
(3*sum [ j 2 , 1 , 3 , sum [ j 1 , 1 , 3 , conj [Yu [j 1 , j 2] ] *Yu [ j 1 , j 2] ] ] *Yu [il , i2] ) /2 + 
(3*sum [ j 2 , 1 , 3 , sum [ j 1 , 1 , 3 , conj [Yu [j 2 , j 1] ] *Yu [ j 2 , j 1] ] ] *Yu [il , i2] ) 11 - 
( 13*gl~2*Yu [il , i2] ) /15-3*g2~2*Yu [il , i2] - (16*g3~2*Yu [il , i2] ) /3 

5. Loop Corrections 

SARAH calculates the analytical expressions for the one-loop corrections to the tadpoles and the self energy 
of all particles. These calculations are performed in DR-scheme and in the 't Hooft gauge. In principle, this 
is a generalization of the calculations for the MSSM presented in [l(|. The command to start the calculation 
is 

CalcLoopCorrections [Eigenstates] ; 

As usual, Eigenstates can in the case of the MSSM either be GaugeES for the gauge eigenstates or EWSB 
for the eigenstates after EWSB. If the vertices for the given set of eigenstates were not calculated before, 
this is done before the calculation of the loop contributions begins. 

Conventions. Using the conventions of [Icj . the results will contain the Passarino Veltman integrals listed 
in |Appcndix D.2| The involved couplings are abbreviated by 

• Cp[pl,p2,p3] and Cp [pi ,p2 ,p3 ,p4] for non-chiral, three- and four-point interactions involving the 
particles pi - p4. 

• Cp[pl,p2,p3] [PL] and Cp[pl,p2,p3] [PR] for chiral, three-point interactions involving the fields pi 
-p3. 

The self-energies can be used for calculating the radiative corrections to masses and mass matrices, respec- 
tively. We have summarized the needed formulas for this purpose in |Appcndix D.23] For calculating the 
loop corrections to a mass matrix, it is convenient to use unrotated, external fields, while the fields in the 
loop are rotated. Therefore, SARAH adds to the symbols of the external particle in the interaction an U for 
'unrotated', e.g. Sd — » USd. The mixing matrix associated to this field in the vertex has to be replaced by 
the identity matrix when calculating the correction to the mass matrix. 
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Results. The results for the loop corrections are saved in two different ways. First as list containing the dif- 
ferent loop contribution for each particle. Every entry includes the following information: internal particles, 
generic type of the diagram, numerical factors coming from symmetry considerations and possible charges 
in the loop. The second output is a sum of all contributions, where the generic results of | Appendix D.2 have 



already been inserted. This can afterwards be written as pdf file using the r^TgX output of SARAH. 
The results for the self-energies are saved in Self EnergylLoopList as list of the contributions and in 
Self EnergylLoopSum written as sums. The last one is a two-dimensional array. The first column gives the 
name of the external particle, the entry in the second column depends on the type of the field. For scalars 
the one-loop self energy II(p 2 ) is given, for fermions the one-loop self-energies for the different polarizations 
(£ L (p 2 ),£ fl 2 ), S s (p 2 )) are written in a three-dimensional array, while for vector bosons the transversal 
part of the self energy H T (p 2 ) is shown. 

Also the results for corrections to the tadpoles are saved twice: in Tadpoles lLoopSum [Eigenstates] ex- 
plicitly written as sum and in Tadpoles lLoopList [Eigenstates] as list of the different contributions. 

Examples. 

1. One-loop tadpoles. The radiative correction of the tadpoles due to a chargino loop is saved in 
TadpoleslLoopList [EWSB] [[1]] ; 

and reads 

{bar[Cha] ,Cp [Uhh[{g01}] ,bar [Cha[{gll>] ] ,Cha[{gIl}]] ,FFS, 1,1/2} 

The meaning of the different entries is: (i) a chargino (Cha) is in the loop, (ii) the vertex with an 
external, unrotated Higgs (Uhh) with generation index gOl and two charginos with index gll is needed, 
(iii) the generic type of the diagram is FFS, (iv) the charge factor is 1, (v) the diagram is weighted by 
a factor ^ with respect to the generic expression (see Appendix D.2 1. 



The corresponding term in Tadpoles lLoopSum [EWSB] is 

4*sum [gll , 1 , 2 , AO [Mass [bar [Cha [{gll}] ] ] ~2] * 

Cp [phid , bar [Cha [{gll}] ] , Cha [{gll}] ] *Mass [Cha [{gl 1}] ] ] 

One-loop self-energies. 

(a) The correction to the down squark matrix due to a four-point interaction with a pseudo scalar 
Higgs is saved in Self EnergylLoopList [EWSB] [ [1 , 12]] and reads 

{Ah,Cp[conj [USd[{g01}]] ,USd[{g02}] ,Ah[{gIl}] ,Ah[{gIl}]] ,SSSS, 1,1/2} 
This has the same meaning as the term 
-sum [gll , 1 , 2 , AO [Mass [Ah [{gl 1}] ] ~2] * 

Cp[conj [USd[{g01}]] ,USd[{g02}] ,Ah[{gIl}] , Ah [{gll}] ] ] /2 
in Self EnergylLoopSum [EWSB] . 

(b) Corrections to the Z boson are saved in Self EnergylLoopList [EWSB] [ [15] ] . An arbitrary entry 
looks like 

{bar[Fd], Fd, Cp[VZ, bar [Fd [{gll}]] , Fd[{gI2}]], FFV, 3, 1/2} 

and corresponds to 

(3*sum[gll, 1, 3, sum[gI2, 1, 3, 

H0[p~2, Mass [bar [Fd [{gll}]]] '2, Mass [Fd [{gI2}] ] ~2] * 
(conj [Cp[VZ,bar[Fd[{gIl}]] ,Fd[{gI2}]] [PL]]* 
Cp[VZ,bar[Fd[{gIl}]] ,Fd[{gI2}]] [PL] + 
conj [Cp[VZ,bar[Fd[{gIl}]] ,Fd[{gI2}]] [PR]]* 
Cp[VZ,bar[Fd[{gIl}]] ,Fd[{gI2}]] [PR]) + 
2*B0 [p~2 ,Mass [bar [Fd [{gll}] ] ] ~2 ,Mass [Fd [{gI2}] ] "2] * 
Mass [bar [Fd [{gll}] ] ] *Mass [Fd [{gI2}] ] * 
Re [Cp [VZ , bar [Fd [{gl 1}] ] , Fd [{gI2}] ] [PL] * 

Cp [VZ , bar [Fd [{gl 1}] ] , Fd [{gI2}] ] [PR] ) ] ] ) /2 
in Self EnergylLoopListSum [EWSB] . 
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6. Definition of models 



All information of the different models is saved in three different hies, which have to be in one directory 
PackageDirectory/Models/ModelName/ 

The three files are: one model file with the same name as the directory (ModelName .m), a file containing 
additional information about the particles of the model (particles .m) and a file containing additional 
information about the parameters of the model (parameters .m). Only the first file is really necessary for 
calculating the Lagrangian and to get a first impression of a model. However, the other two files are needed 
for defining properties of parameters and particles and for producing an appropriate output. 
In addition, it is possible to include LesHouches spectrum files in the model directory [l6|. These can be 
read with the command ReadSpectrum ["Filename"] . If spectrum files are located in another directory, the 
complete path has to be added to the filename to read them. 

6.1. The model file 

The model file contains the following parts: first, the gauge structure and the particle content are given. 
Afterwards, the matter interactions are defined by the superpotential. These are general information which 
must always be apparent. 

This part is followed by the definition of the names for all eigenstates (NameOf States). For these eigen- 
states, several properties can be defined using the corresponding DEFINITION statement: decomposition 
of complex scalars in scalar, pseudo scalar and VEV (DEFINITION [States] [VEVs]), rotations in the mat- 
ter (DEFINITION [States] [MatterSector] ) and gauge sector (DEFINITION [States] [GaugeSector] ), the 
corresponding gauge fixing terms (DEFINITION [States] [GaugeFixing] ), the flavor decomposition of fields 
(DEFINITION [States] [Flavors]) and possible phases of fields (DEFINITION [States] [Phases]). New cou- 
plings can be added and existing couplings can be changed by hand (DEFINITION [States] [Additional]). 
Afterwards, optionally the particles are states, which should be integrated out or deleted. At the end, the 
Dirac spinors have to be built out of Weyl spinors, a spectrum file can be defined and a choice for an 
automatic output can be made. 

6.1.1. Vector and chiral superfields 

SARAH supports all SU(N) gauge groups. The gauge sector in SARAH is defined by adding a vector 
superfield for each gauge group to the list Gauge in the model file, e.g. 

Gauge[[3]]={G, SU[3], color, g3, False}; 

The different parts define the names of the superfield, of the gauge group and of the gauge coupling. In 
addition, the dimension of the gauge group is given. The last entry states, if the gauge indices should be 
implicit or explicit. The name of the gaugino component of the vector superfield starts with f and the vector 
boson with V, i.e. based on the above definition, the gluino is called f G and the gluon VG. The ghost field is 
named using g, i.e. gG. 

Chiral superfields are defined by using the list Fields, e.g. 
Fields [[1]] = {{uL, dL}, 3, q, 1/6, 2, 3}; 

Fields [[5]] = {conj [dR] , 3, d, 1/3, 1, -3}; 

The first entry defines the names used for the component fields, then the number of generation and the 
name for the superfield follows. The automatically created name of the fermionic component starts with F 
and the scalar one with S: the squarks are named SuL, SdL or SdR, while the quarks are FuL, FdL and FdR. 
After the superfield name, the representation with respect to the gauge groups defined by Gauge is assigned. 
The transformation of an irreducible representation r under a given gauge group is in most cases fixed by 
its dimension D. Therefore, it is sufficient to assign only D if it is unique. Otherwise, the Dynkin labels of 
r have to be given as additional input. 
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For all gauge groups, the generators of all appearing representations are needed in order to write the kinetic 
part of the Lagrangian and the D-terms. In this context, all generators for non-fundamental representations 
are written as tensor products in SARAH. Furthermore, the eigenvalues C2(r) of the quadratic Casimir as well 
as the Dynkin index I(r) are needed. The performed steps to obtain these results are given in | Appendix C[ 

6.1.2. Superpotential 

The superpotential is defined in a compact form using the variable SuperPotential: 

SuperPotential = {{{Coefficient, Parameter, (Contraction)}, 

{Particle 1, Particle 2, Particle 3} >, ...} 

Each term of the superpotential is defined by two lists. The second list assigns all involved fields by using the 
superfield names. The first list is two- or three-dimensional. It defines a numerical coefficient and the name 
of the coupling. The gauge and generation indices of the involved supcrficlds arc automatically contracted 
by SARAH. The used contraction can be displayed via 

ShowSuperpotentialContractions ; 

Sometimes, there are more possibilities to contract all indices. Thus, it is possible to fix the contraction of 
each term using the third entry of the first list. 

6.1.3. Symmetry breaking and rotations 

Rotations for all matter and gauge fields as well as the decomposition of scalar fields into their scalar 
component, pseudo scalar components and VEVs can be performed. All appearing coefficients as well as the 
names of the rotation matrices to parametrize this change of the basis can be chosen individually. Besides, 
it is possible to decompose a field carrying a generation index into its different flavors in order to treat them 
separately. Afterwards, the complete Lagrangian for the new set of eigenstates is calculated. 
We give here again only the generic syntax for the different tasks and refer to 16.1.51 for a discussion of the 
MSSM. The definition of rotations in the matter sector has either the form 

{Old Basis, {New Basis, Mixing Matrix}} 

or 

{{l.Qld Basis, 2. Old Basis}, 

{{l.New Basis, 1. Matrix}, {2. New basis ,2 .Matrix}}} 

depending on whether the mass matrix is hermitian or not. The decomposition of scalar fields is done via 

{Scalar, {{VEV, l.Coeff.}, {Pseudoscalar , 2 . Coef f . } , {Scalar ,3 . Coef f . }} ; 

Finally, the syntax for the flavor decomposition of fields is 

{Field, {Name for 1 . Generation, Name for 2. Generation, ...}}; 

In principle, these steps can be repeated as often as needed. Therefore, it is no problem to go first to the 
SCKM basis and afterwards to the mass eigenbasis. GUT theories involving several symmetry breakings 
can be treated in the same way. The information of all intermediate steps is saved. Hence, it is possible 
to calculate the vertices or masses of all eigenstates without the necessity of a new model file or a new 
evaluation of the model. 
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6.1. 4-- Effective or non- super symmetric theories 

It is easy in SARAH to integrate particles out of the spectrum to get an effective theory, or just to delete 
some particles to get a non-supersymmetric limit of the model. Integrating out particles happens by adding 
them to the list IntegrateOut 

IntegrateOut = {Particle 1, Particle 2, . ..} 

If particles are integrated out, all higher dimensional operators up to dimension 6 are calculated. Delet- 
ing particles can be done in the same way as integrating them out. The corresponding list is called 
DeleteParticle. The difference is that there are no effective operators calculated. Deleting is therefore, of 
course, faster and should be used if the higher dimensional operators are not needed. 

6.1.5. Example: Model file for the MSSM 

We show in the following the implementation of the MSSM in SARAH. Our conventions are discussed in 
|Appcndix E| 

1. The gauge sector is U(l) x SU(2) x SU(3) and is just defined by adding the corresponding vector 
superfields. 

Gauge [[1]]={B, U[l] , hypercharge, gl, False}; 
Gauge [ [2] ] ={WB , SU[2], left, g2, True}; 

Gauge[[3]]={G, SU[3], color, g3, False}; 

2. The doublet superfields q, I, H d , H u are added by 

Fields [[1]] = {{uL, dL} , 3, q, 1/6, 2, 3}; 

Fields[[2]] = {{vL, eL} , 3, 1, -1/2, 2, 1}; 

Fields [[3]] = {{HdO, Hdm} , 1, Hd, -1/2, 2, 1}; 

Fields [[4]] = {{Hup, HuO} , 1, Hu, 1/2, 2, 1}; 

3. The singlet superfields d, u, e are added by 

Fields [[5]] = {conj [dR] , 3, d, 1/3, 1, -3}; 
Fields[[6]] = {conj [uR] , 3, u, -2/3, 1, -3}; 
Fields [[7]] = {conj [eR] , 3, e, 1, 1, 1}; 

4. The superpotential of the MSSM is 

W = qY u uH u - qY d dH d - lY e eH d + fiH u H d (3) 

and given in SARAH by 

SuperPotential = { {{1, Yu},{u,q,Hu}}, {{-l,Yd},{d,q,Hd}}, 

{{-l,Ye},{e,l,Hd}}, {{1 , \ [Mu] } , {Hu.Hd}} }; 

5. There are two different sets of eigenstates: the gauge eigenstates before EWSB and the mass eigenstates 
after EWSB. The internal names are 

NameOf States={GaugeES, EWSB}; 

6. The gauge fixing terms for the unbroken gauge groups are 

DEFINITION [GaugeES] [GaugeFixing] = 
{ {Der[VWB], -1/(2 RXi[W])}, 
{Der[VG] , -1/(2 RXi [G] ) }}; 

This has the same meaning as 

£gf = —^W^ - ±\ d M g^\ 2 (4) 
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7. The vector bosons and gauginos rotate after EWSB as follows 
DEFINITION [EWSB] [GaugeSector] = 

{ {VWB, {l,{VWm, 1/Sqrt[2]>, {conj [VWm] , 1/Sqrt [2] }} , 

{2,{VWm,-I/Sqrt[2]>, {conj [VWm] , I/Sqrt [2] }} , 

{3 , {VP , Sin [ThetaW] } , {VZ, Cos [ThetaW] }» , 
{VB, {1, {VP, Cos [ThetaW]}, {VZ, -Sin [ThetaW] »> , 
{fWB, {l,{fWm, 1/Sqrt[2]}, {fWp 1/Sqrt[2]», 

{2,{fWm,-I/Sqrt[2]>, {fWp, I/Sqrt[2]}}, 

{3,{fW0, 1}}} 

}; 

This is the common mixing of vector bosons and gauginos after EWSB, sec | Appendix E3| 

8. The neutral components of the scalar Higgs receive VEVs Vd/v u and split in scalar and pseudo scalar 
components according to |Appendix E.4 This is added to SARAH by 



DEFINITION [EWSB] [VEVs] = 

{{SHdO, {vd, 1/Sqrt[2]>, {sigmad, I/Sqrt [2] } , {phid, 1/Sqrt[2]}}, 
{SHuO, {vu, 1/Sqrt[2]>, {sigmau, I/Sqrt [2] }, {phiu, 1/Sqrt[2]»>; 

9. The particles mix after EWSB to new mass eigenstates 
DEFINITION [EWSB] [MatterSector] = 
{{{SdL, SdR }, {Sd, ZD}}, 

{{SuL, SuR }, {Su, ZU}}, 

{{SeL, SeR }, {Se, ZE}}, 

{{SvL }, {Sv, ZV}}, 

{{phid, phiu }, {hh, ZH}}, 

{{sigmad, sigmau }, {Ah, ZA}}, 
{{SHdm, conj[SHup] }, {Hpm.ZP}}, 
{{fB, fWO, FHdO, FHuO}, {LO, ZN}}, 
{{{fWm, FHdm}, {fWp, FHup}}, {{Lm.Um}, {Lp.Up}}}, 
{{{FeL}, {conj [FeR] }} , {{FEL , ZEL} , {FER , ZER}}} , 

{{{FdL}, {conj [FdR]}}, {{FDL , ZDL} , {FDR , ZDR}}} , 

{{{FuL}, {conj [FuR]}}, {{FUL,ZUL},{FUR,ZUR}}} }; 

This defines the mixings to the mass eigenstates described in |Appcndix E.6| 

10. The new gauge fixing terms after EWSB are 
DEFINITION [EWSB] [GaugeFixing] = 

{{Der[VP] , - 1/(2 RXi[P])}, 

{Der[VWm]+ I Mass [VWm] RXi [W] Hpm[{l}], - l/(RXi[W])}, 

{Der[VZ] - Mass[VZ] RXi [Z] Ah[{l}], - 1/(2 RXi[Z])}, 

{Der[VG] , - 1/(2 RXi [G] ) }} ; 

Based on this definition, A® and Hf- are recognized in all calculations as Goldstone bosons. 

11. No particles should be integrated out or deleted 
IntegrateOut={} ; 
DeleteParticles={} ; 

12. The Dirac spinors for the mass eigenstates are 

dirac[[l]] = {Fd, FdL, FdR}; 

dirac [[2]] = {Fe, FeL, FeR}; 

dirac [[3]] = {Fu, FuL, FuR}; 

dirac [[4]] = {Fv, FvL, 0}; 

dirac [[5]] = {Chi, LO, conj[L0]}; 



dirac [[6]] = {Cha, Lm, conj [Lp] } 
dirac [[7]] = {Glu, fG, conj [f G] } 
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6.2. Parameter and particle files 

6.2.1. Parameter file 

Additional properties and information about the parameters and particles of a model are saved in the 
files parameters .m and particles. m. An entry in the parameter file looks like 

{Yu, { LaTeX -> "Y~u", 
Real -> True, 
Form -> Diagonal, 
Dependence -> None, 
Value -> None, 
LesHouches -> Yu 
}} 

and contains information about the numerical value (Value — > number), the position in a LesHouches 
accord file (LesHouches — »■ position) or the dependence on other parameters (Dependence — > equation). 
Also simplifying assumptions can be made: it can be defined that parameters contain only real entries (Real 
— >■ True) or that the parameter is diagonal (Form — > Diagonal). Also a F/Tj^X name can be given (LaTeX — > 
name). Furthermore, the GUT normalization can be assigned (GUTnormalization — > value) for the gauge 
couplings of an U(l ) gauge groups. 

6.2.2. Particle file 

The particles file contains entries like 

{Su , { RParity -> -1, 

PDG -> {1000002,2000002,1000004,2000004,1000006,2000006}, 

Width -> Automatic, 

Mass -> Automatic, 

FeynArtsNr -> 13, 

LaTeX -> "\\tilde{u}" , 

DutputName -> "urn" }}, 

and defines properties of all particles such as the i?-parity (RParity — > number) or the mass (Mass — > value 
or Automatic). Automatic means that for the output for FeynArts or CalcHep not a fixed numerical value 
is used, but that the masses are calculated using tree-level relations. In addition, the PDG code is given 
(PDG — > number), the number for the particle class used in the FeynArts model file can be fixed (FeynArts 
— > number) and the name in FaTjtX form is given (LaTeX — > name). If a CalcHep or CompHep model file 
should be written, it is also helpful to define an appropriate name in this context (DutputName — > name). 

6.2.3. Global definitions 

It is also possible to define global properties for parameters or particles which are present in more than 
one model file. These properties are afterwards used for all models. The global definitions are saved in the 
files particles. m and parametes.m directly in the main model directory. For each parameter or particle, 
an entry like 

{{ Descriptions -> "Up Squark" , 

RParity -> -1, 

PDG -> {1000002,2000002,1000004,2000004,1000006,2000006}, 

Width -> Automatic, 

Mass -> Automatic, 

FeynArtsNr -> 13, 

LaTeX -> "\\tilde{u}" , 

DutputName -> "urn" }}, 
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can be added. In particular, the entry Description is important. This should be an unique identifier for 
each particle or parameter. This identifier can later on be used in the different files of the different models, 
e.g. 

{Su , { Descriptions -> "Up Squark"}}, 

Of course, it is also possible to overwrite some of the global definitions by defining them locally, too. For 
instance, to use u instead of um as output name in a specific model, the entry should be changed to 

{Su , { Descriptions -> "Up Squark", 
OutputName -> "u" }}, 

in the corresponding particle file of the model. 



7. Verification 

Tree-level results. We have checked the model files generated with SARAH for the MSSM against the existing 
files of FeynArts and CalcHep. The checks happened as well at vertex level as for complete processes. 
We have compared the numerical value of each vertex for different sets of parameters and all possible 
combinations of generations (more than 5000). In addition we have calculated several 1 — > 2 and 2^2 
processes with the old and new model files. Similar checks have been done for the vertices of the NMSSM 
against the model file of CalcHep. More information about the verification of the tree-level results is given 



11] 



One-loop self-energies and tadpole equations. We have compared the analytical expressions of the self- 
energies calculated by SARAH for the MSSM with the results of [ltj. In addition, we have compared the 
numerical values for the one-loop corrected masses with the results of SPheno [2(| . Furthermore, the results 
for the NMSSM were numerically checked against the routines provided by the authors of 21]. Both results 
were in complete agreement. A detailed discussion about the results for the NMSSM obtained by SARAH is 



also given in [22 



Verification of the calculated RGEs. We have compared the analytical results for the one- and two-loop 
RGEs calculated by SARAH for the MSSM with [§] and for the NMSSM with 0. The only difference has 
been in the NMSSM the two- loop RGE of A\. A second calculation by authors of [23| have confirmed the 
first result of SARAH. 

In addition, we have checked a model containing non-fundamental representations: the SU (5) inspired 
Seesaw II model of [24| and [25[ . It is known that there are discrepancies of the RGEs given in these two 



papers. Our result fully agrees with [24 



Also numerically checks have been done by comparing the RGEs for the MSSM with the RGEs implemented 
in SPheno 20]. Both sets of RGEs are in full agreement. 
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Appendix A. Gauge anomalies 



Before SARAH starts the calculation of the Lagrangian it checks the model for the different triangle 
anomalies. These anomalies can involve diagrams with three external gauge bosons belonging to the same 
U(l) or SU(N) gauge group. To be anomaly free, the sum over all internal fermions has to vanish 

U(l)l : £lf = 0, (A.l) 

n 

SU(N)I : ^2 Tl ( T n T n T n) = ■ (A.2) 
n 

We label the different gauge groups with the indices i,j, k. Y£ is the charge of particle n under the abelian 
gauge group i while is the generator with respect to a non-abelian gauge group. 

Combinations of two different gauge groups are possible, if one group is an U(l). Hence, another condition 
for the absence of anomalies is 

U(l)i x SU(N)* : J2 Y n Tr(T^T^) = . (A.3) 

n 

If more than one U(l) gauge group are present, anomalies can be generated by two or three different U(l) 
gauge bosons as external fields, too. Therefore, it has to be checked, that 

U(l)i x [/(l), 2 : J2 Y n Y f = °> ( A -4) 

n 

U(l) t xU(l) 3 xU(l) k : ^y^Y n fe =0 (A.5) 

n 

holds. In addition, it has to be checked that there is an even number of SU(2) doublets. This is necessary 



for a model in order to be free of the Witten anomaly [26j. If one condition is not fulfilled, a warning is 
given by SARAH but the model can be evaluated anyway. 

Appendix B. Calculation of the Lagrangian of supersymmetric models 

We describe in this section the calculation of the complete Lagrangian for a supersymmetric model based 
on the superpotential and the gauge structure. 

Interactions of chiral superfields. If we call the superpotential for a given theory W and use (f>i for the scalar 
and ipi for the fermionic component of a chiral supermultiplet, the matter interactions can by derived by 

£ Y = --W ij ip i %l> j +b.c. , £ F = F^Fi + h.c. (B.l) 

with 

S 2 SW 

W 13 = — — W and F l = -W* 1 = — . (B.2) 

ocpiocpj bepi 



The first term of eq. (|B.1[) describes the interaction of two fermions with one scalar, while the second term 
forms the so called F-terms which describe four-scalar interactions. 

Interactions of vector superfields. We name the spin-^ component of a vector supermultiplet A and the 
spin-1 component A^. The most general Lagrangian only involving these fields is 

L = ~F^ v F' xva - iX^a-^D^X" (B.3) 



with the field strength 



F% = d,Al - d v Al + gf abc A\Al , (B.4) 
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and the covariant derivative 

D^X a = d^X a +gf abc Al\ c . (B.5) 

Here, f abc is the structure constant of the gauge group. Plugging eq. (|B.4|) in the first term of eq. (|B.3I) 
leads to self-interactions of three and four gauge bosons 

L V = -\{d^A a v - d v Al)gf abc A^ b A^ - ^g 2 {f abc A b ll Al){f ade A^ e A v ' e ) . (B.6) 



The second term of eq. (|B.3[) describes the interactions between vector bosons and gauginos. 

Super symmetric gauge interactions. The parts of the Lagrangian with both chiral and vector superfields are 
the kinetic terms for the fermions and scalars 

<C«„ = -D^^D^ - i^a^D^i (B.7) 

as well as the interaction between a gaugino and a matter fermion and scalar 

L GFS = -V2. 9 (0*T<V)A a + h.c. . (B.8) 

Here, T a are the fundamental generators of the gauge group. Furthermore, the covariant derivatives are 

D^i = d^i-igA^T^i, (B.9) 
= d^ + igA^T-y , (B.10) 
D^i = d^ipi — igA a ^(T a ^)i , (B.ll) 

In addition, the D-terms are defined by 

L D = ^D a D a . (B.12) 

The solution of the equations of motion for the auxiliary fields leads to 

D a = - 5 (0*T» . (B.13) 

Soft-breaking terms. SUSY must be a broken. This can be parametrized by adding soft-breaking terms to 
the Lagrangian. The possible terms are the mass terms for all scalar matter fields and gauginos 

&SB = -ml i <l>i<t>;-±Mx i XiXi (B.14) 

as well as soft-breaking interaction corresponding to the superpotential terms 

&Soft,w = T(j)i(j)j4> k + Bfafa + S4>i . (B.15) 

Gauge fixing terms and ghost interactions. The Lagrangian of a theory without further restrictions is in- 
variant under a general gauge transformation. This invariance leads to severe problems in the quantization 
of the theory as can be seen in the divergence of functional integrals. Therefore, it is in necessary to add 
gauge fixing terms to break this gauge invariance. 
The general form of the gauge fixing Lagrangian is 



GF 



^Ei/^i 2 - ( R16 ) 



f a can be a function of partial derivatives of a gauge boson and a Goldstone boson. The corresponding 
ghost terms of the ghost fields fj and r\ are 

^Ghost = -fja{Sf a ) ■ (B.17) 
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Here, S assigns the operator for a BRST transformation. For an unbroken gauge symmetry, the gauge fixing 
terms in the often chosen -gauge are 



£e F = - w Y,( d "K) 2 ■ (B.18) 

Here, are the gauge boson of the unbroken gauge group. It is often common to choose a distinct value for 
R^. The most popular gauges are the unitary gauge R^ —¥ oo and the Feynman-'t Hooft-gauge R^ = 1. For 
broken symmetries, the gauge fixings terms are chosen in a way that the mixing terms between vector bosons 
and scalars disappear from the Lagrangian. Therefore, the common choice for the gauge fixing Lagrangian 
for theories with the standard model gauge sector after EWSB is 

Zgfm, = (d^f - (d»Z, + izM z G a ) 2 + (&>W+ + tw+M w G+) 2 . (B.19) 

Here, G° and G + are the Goldstone bosons, which build the longitudinal component of the massive vector 
bosons. 



Appendix C. Calculation of Group Factors 

SARAH supports not only chiral superficlds in the fundamental representation but in any irreducible 
representation of SU(N). In most cases, it is possible to fix the transformation properties of the chiral 
superfield by stating the corresponding dimension D. If the dimension is not unique, also the Dynkin labels 
are needed. For calculating kinetic terms and D-terms, it is necessary to derive from the representation 
the corresponding generators. Furthermore, the eigenvalues C2 of the quadratic Casimir for any irreducible 
representation r 

T a T a 4>(r)=C 2 {r)4>(r) (C.l) 

as well as the Dynkin index I 

Tr(T a T b )4>(r) = I5 ab <j>{r) (C.2) 

are needed for the calculation of the RGEs. All of that is derived by SARAH due to the technique of Young 
tableaux. The following steps are evolved: 

1. The corresponding Young tableau fitting to the dimension D is calculated using the hook formula: 

D = TU^, (C3) 
di is the distance of the i. box to the left upper corner and hi is the hook of that box. 

2. The vector for the highest weight A in Dynkin basis is extracted from the tableau. 

3. The fundamental weights for the given gauge group are calculated. 

4. The value of C2M is calculated using the Weyl formula 

C 2 (r) = (A,A + p), (C.4) 

p is the Weyl vector. 

5. The Dynkin index I(r) is calculated from C2(r). For this step, the value for the fundamental repre- 



sentation is normalized to |. 



J(r)=C 2 (r)^, (G5) 



with D(G) as dimension of the adjoint representation. 
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6. The number of co- and contra- variant indices is extracted from the Young tableau. With this infor- 
mation, the generators are written as tensor product. 

The user can calculate this information independently from the model using the new command 
ChecklrrepSUN [Dim.N] 

Dim is the dimension of the irreducible representation and N is the dimension of the SU (N) gauge group. 
The result is a vector containing the following information: (i) repeating the dimension of the field, (ii) 
number of covariant indices, (iii) number of contravariant indices, (iv) value of the quadratic Casimir C2(r), 
(v) value of the Dynkin index I(r), (vi) Dynkin labels for the highest weight. 

Examples. 

1. Fundamental representation. The properties of a particle, transforming under the fundamental 
representation of SU(3) are calculated via ChecklrrepSUN [3,3] . The output is the well known result 

{3, 1, 0, 4/3, 1/2, {1, 0}} 

2. Adjoint representation. The properties of a field transforming as 24 of SU(5) are calculated by 
ChecklrrepSUN [24, 5] . The output will be 

{24, 1, 1, 5, 5, {1, 0, 0, 1» 

3. Different representations with the same dimension. The 70 under SU(5) is not unique. There- 
fore, ChecklrrepSUN [{70 , {0, 0, 0, 4}}, 5] returns 

{70, 0, 4, 72/5, 42, {0, 0, 0, 4» 

while ChecklrrepSUN [{70, {2, 0, 0, 1}}, 5] leads to 

{70, 2, 1, 42/5, 49/2, {2, 0, 0, 1}} 



Appendix D. Conventions and generic expressions 

Appendix D.l. Renormalization group equations 

We summarize in this section the used equations for the calculation of the one- and two-loop RGEs in 
SARAH. These equations are extensively discussed in literature, see e.g. 
For a general N = 1 supcrsymmctric gauge theory with superpotential 

W((f>) = Li4>i + ^n iS Mi + ^V !k o,o,u, , (D.l) 



the soft SUSY-breaking scalar terms are given by 

^soft = (s 1 ^ + + \h ijk <i>i<t>i<t>k + ccj + (m 2 ) + ~M\A a A a 

The anomalous dimensions are given by 



(D.2) 



ir=-Y ipq Y^-2Sig 2 C(i), (D.3) 

7f W = - \yr mn Y npq Y pqr Y m ^ + g 2 Y tpq Y^[2C(p) - C{i)\ 

+ 25lg i [C{i)S{R) + 2C(i) 2 - 3C(G)C(i)} , (D.4) 

and the /3-functions for the gauge couplings are given by 

/3« - 9 3 [S(R) - 3C(G)] , (D.5) 

Pf ] =9 5 |-6[C(G)] 2 + 2C(G)S(R) + AS(R)C(R)} - g 3 Y^ k Y l]k C(k)/D(G) . (D.6) 
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Here, C(i) is the quadratic Casimir for a specific superfield and C(R),C(G) are the quadratic Casimirs for 
the matter and adjoint representations, respectively. D(G) is the dimension of the adjoint representation. 
The corresponding RGEs are defined as 



d 

di 9 



1 



16tt 2 



(16tt 



2\2 * 9 



(D.7) 



Here, we used t = InQ, where Q is the renormalization scale. The /3-functions for the superpotential 
parameters can be obtained by using superfield technique. The obtained expressions are 



Y 

H ip 
L p 



i j i' 



1 



(i)fc_ 



1 



16tt 2 lp 
1 



16^ 2 ,p 

1 



(2U- 



7i iW + 



(16^ 2 ) 2 lp 
1 



7«' + 



(16tt 2 ) 2 
1 



(2)j 



+ (fc o z) + (fc <-> j) 

+ (j *) , 



16tt 2 lp ' (16tt 2 ) 27p 
The expressions for trilinear, soft-breaking terms are 



(2)i 



16^ 2 



(i) 



ijk 



(16^ 2 ) 2 



(2) 



with 



(1) 



ijk X 



+ (k •H- i) + (A; j) , 



ft 



(2) 



ijk 



1 



hWV, v n Plv V 

11 1 Iran 1 1 pqr 1 



rark 



rmrk 



_ V l 3 l V, V n P"V umrk _ yijl v unpqy „ 
1 1 Iran- 1 1 pqr 11 1 1 twin 11 1 pqr 1 

+ (h ijl Yi pq Y pqk + 2Y ijl Yi pq h pqk - 2MY^ l Yi pq Y pqk ) g 2 [2C(p) - C(k)\ 
+ (2h ijk - 8MY ijk ) g A [C(k)S(R) + 2C(k) 2 - 3C(G)C(k)] 
+ (k <h> i) + (k «• j) . 



For the bilinear soft-breaking parameters, the expressions read 



dt 16tt 2 



+ 



1 



(16^ 2 ) 2 



?(2) 



with 



/5 



A 



(2) 



— h il Y, vranj , }_yiflv 
"2 Inin 1 ~ 2 1 Iran 

+ {l<r> j) , 



+ ^^ m nft m ^ - 2 (&« - 2M/i«) 5 2 C7(i) 



Lmry ypqn 
mnU 1 pq r 1 



— ZVWV, ,,' mr Y hPQ n — n il Vi h n ^Y Y' mr '3 
2 J 1 ImnH' ±pq r ii> p J-lmn 11, 1 pqr 1 

- ii il Y lmn Y npq Y pqr h mrj + 2Y^ l Y lpq (b pq - l i pq M)g 2 C{p) 

+ (b a Y lpq Y p <» + 2^ l Y lpq hP qi - 2^ l Y lpq Y^M) g 2 [2C(p) - C(i)] 
+ (2b^ - 8/i y M) g A [C(i)S(R) + 2C(i) 2 - 3C(G)C(i)] 
+ (» «• i) , 
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(D.8) 
(D.9) 
(D.10) 

(D.ll) 



(D.12) 



(D.13) 



(D.14) 



(D.15) 



(D.16) 



Finally, the RGEs for the linear soft-breaking parameters are 



d , 
dt 



1 



16tt 2 



+ 



(16tt 



2\2 



,(2) 



with 



4 1} 



4 2) 



— -Y 1 n Y. p i n S p + L p Y p i n h 



tin 



\1 Y kin & 



2Y lKp {m% m + h lkl B. 



kl , 



(D.17) 



(D.18) 



=2g 2 C(l)Y lkl Y pkl S p - -Y^Y qst Y lst Y pkl S p - 4g 2 C(l)(Y m M - h lkl )Y pkl L p 

- [Y^Y qst h lst Y pM + h^Y qst Y lst Y pkl ]L p - 4g 2 C(l)Y jnl (fx nl M - B nl )^ 

- [Y rnq h« st Y lstf i nl + Y jnq Y« st Y lst B nl ]^ + Ag 2 C{l){2Y lkl m \M\ 2 - Y m B kl M 

- h ikl m M* + h m B kl + Y ipl (m 2 ) k m + Y* kp (m 2 ) l pm ) 

- Y^Y qst h lst B kl + h^Y qst Y lst B kl + Y^h qst h lst m + h ik "h qst Y lst m 
+ Y^(m 2 ) k Y qst Y lst m + Y^Y qst Y pst (m 2 ) l pm + Y lkp (m 2 );Y qst Y lst m 
+ 2Y^Y qsp (m 2 ) p Y lst m . 



(D.19) 



With these results, the list of the /3-functions for all couplings is complete. Now, we turn to the RGEs for 
the gaugino masses, squared masses of scalars and vacuum expectation values. The result for the gaugino 
masses is 



dt 



1 



16tt 2 



Pm 



1 



(16tt 2 ) : 



— B {2) 



with 



Pm 

o(2) 

Pm 



--g 2 [2S(R) - 6C(G)] M , 



= 5 4 |-24[C(G)] 2 + 8C(G)S(R) + 16S(R)C(R)} M 

+ 2g 2 [h^ k - MY i3k ] Y ljk C(k)/D(G) 

The one- and two-loop RGEs for the scalar mass parameters read 



d 

dt m ^ 



16tt 2 



(i) 



(16^ 2 ) 2 



(2) 



(D.20) 

(D.21) 
(D.22) 

(D.23) 
(D.24) 
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with 



,(2) 



\ =\Y lpq Y^(mX + l -Y^Y pqn {m 2 ) n t + 2Y lpq Y^ { m 2 f r 
+ h lpq b? pq - 85jMM^g 2 C(i) + 2g 2 t^Tr[t A m 2 } , 



(D.25) 



1, O.J, 



— —(m' z XY, Y mr 3V YP1 n _ —(m 2 \ 3 Y lmn Y YP<l r Y 
2^ Ji^lmn- 1 1 pqr ± 9 ± mri ± 1 pqn 

- Y llm Y^ m (m 2 ) l r Y npq Y^ - Y llm Y^ m (m 2 ) r n Y rpq Y l P" 

- Y Um Y^(m 2 ) l n Y pqr Y^ m - TY llm Y^Y npq Y m ^{m 2 ) q r 
_ Y-, vi ln h h m Pi — h , hi ln Y Y m P1 

± %lm ± ll npq 11 •Him" j 1 npq 1 

— h-, V3 ln V h' m P ( l — Y, h3 ln h V m P1 
"•ilm 1 ± npq llj J -iLm li '^npq 1 

I 



{m 2 )\Y lpq Y^ + Y lpq Y lpq (m 2 )\ + 4Y lpq Y^ 1 (m 2 )* + 2h ipq V™ 

g 2 [C(p) + C(q) - C(i)} 
2g 2 tf\t A m 2 ) l r Y lpq Y rpq + 8 5 4 tf j Tv[t A C(r)m 2 ] 



2h ipq Y jpq M - 2Y ipq h jpq M^ + AY lpq Y jpq MM^ 



+ 8{g 4 MM^ [2AC(i)S{R) + 48C(i) 2 - 72C{G)C{i)} 
+ 8Sjg i C(i)(TT[S(r)m 2 } - C{G)MM' i ) . 



(D.26) 



The RGEs for a VEV v l are proportional to the anomalous dimension of the chiral superfield whose scalar 
component receives the VEV 



dt 



1 



16tt 



7 (1)i + 

2 IP ~ 



1 



(167T 2 ) 2 



(2)» 



(D.27) 



Appendix D.2. One-loop amplitudes for one- and two-point functions 

We used for the calculation of the one-loop self-energies and the one-loop corrections to the tadpoles in 
DR-scheme the scalar functions defined in [lOj. The basis integrals are 



Ao(m) = 16tt 2 Q 4 -" / , ^ 9 1 
J i (2n) n q z — m 

B (p,m l7 m 2 ) = 16tt 2 Q 4 -" / 

J «(2tt) 



m 2 + ie 



1 



q 2 — m\ + ie 



(q — p) 2 — m 2 , + ie 



(D.28) 
(D.29) 



with the renormalization scale Q. The integrals are regularized by integrating in n = 4 — 2e dimensions. 
The result for Aq is 



. / s 9/1 ,771 

A (m) = m 2 I -, + 1-ln^ 
where 1/e = 1/e — 7_e + ln47r. The function _Bq has the analytic expression 



1 



B (p,m 1 ,m 2 ) = - - In ( — ) - f B (x+) - /b(x-) , 



(D.30) 



(D.31) 



with 



x± 



s ± \/s 2 — 4p 2 (m 2 — ie) 
V 



/ B (x) = ln(l-x)-a;ln(l-a;" 1 )-l , (D.32) 



and s = p 2 — m 2 + m\. All the other, necessary functions can be expressed by Aq and Bq. For instance, 



B 1 (p,m ll m 2 ) = 



A Q (m 2 ) - Ao(mi) + (p 2 + m\ - mj)B (p, TOi,m 2 ) 
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(D.33) 



and 



B 2 2(p,m 1 ,m 2 ) = - 




(D.34) 



Furthermore, for the self-energies of vector bosons, it is useful to define 



B 22 (p, mi,m 2 ) 



F (p,m 1 ,m 2 ) 
G (p,m l7 m 2 ) 
H (p,m 1} m 2 ) 



A {mi) - 2A (m 2 ) - (2p 2 + 2m\ - m\)B {p,mi,m 2 ) , 
(p 2 - ml- ml)B a (p,mi,m 2 ) - A (mi) - A (m 2 ) , 
AB 22 (p,mi,m 2 ) + G(p,m 1 ,m 2 ) , 

■B 22 (p,mi,m 2 ) - jAo(mi) - ^-A (m 2 ) 



(D.35) 
(D.36) 
(D.37) 



(D.38) 



In addition, several coefficients are involved: 

• cs is the symmetry factor: if the particles in the loop are indistinguishable, the weight of the contri- 
bution is only half of the weight in the case of distinguishable particles. If two different charge flows 
are possible in the loop, the weight of the diagram is doubled. 

• cc is a charge factor: for corrections due to vector bosons in the adjoint representation this is the 
Casimir of the corresponding group. For corrections due to matter fields this can be, for instance, a 
color factor for quarks/squarks. For corrections of vector bosons in the adjoint representation this is 
normally the Dynkin index of the gauge group. 

• cr is 2 for real fields and Majorana fermions in the loop and 1 otherwise. 

We use in the following T for non-chiral interactions and T^/Tr for chiral interactions. If two vertices are 
involved, the interaction of the incoming particle has an upper index 1 and for the outgoing field an upper 
index 2 is used. 

Appendix D.2.1. One-loop tadpoles 

1. Fermion loop (generic name in SARAH: FFS): 



T = %c s c c m F TAQ(m 2 F ) 



(D.39) 



2. Scalar loop (generic name in SARAH: SSS): 



T = -2c s c c TA Q (m 2 s ) 



(D.40) 



3. Vector boson loop (generic name in SARAH: SVV): 



T = 6c s c c TA (m 2 v ) 



(D.41) 



Appendix D.2.2. One-loop self-energies 
Corrections to fermion. 
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1. Fermion-scalar loop (generic name in SARAH: FFS): 

S s (p 2 ) = m F csc c c R T R T 2 £*B (p 2 ,m 2 F ,m 2 s ) 
E*(p 2 ) = -c s c c c R ^T 1 R T 2 I i*B 1 (p 2 ,m F ,m 2 s ) 

Z L (p 2 ) = -csc c c R ^r]TY Bi(p 2 ,m|,, ml) 

2. Fermion- vector boson loop (generic name in SARAH: FFV): 

S 5 (p 2 ) = -±csc c c R m F TlT 2 £B (p 2 ,m 2 F ,m 2 s ) 
E'V) - -cscccnTlr^B^^j^m 2 ) 
£ L (p 2 ) = -c s c c c R r R T 2 I i*B 1 (p 2 ,rn 2 F ,m 2 s ) 

Corrections to scalar. 

1. Fermion loop (generic name in SARAH: FFS): 

U(p 2 ) = c s c c c R ((Iir 2 '* + r^Go^ 2 , m 2 F , m|) + (rir|* + r^r 2 '*)S (p 2 , m 2 Fl m 2 )) (D.42) 

2. Scalar loop (two 3-point interactions, generic name in SARAH: SSS): 

IV) = csccCRT^Boip 2 ,™, 2 ,,™, 2 ) (D.43) 

3. Scalar loop (4-point interaction, generic name in SARAH: SSSS): 

n(p 2 ) - -c s c c TA (m 2 ) (D.44) 

4. Vector boson-scalar loop (generic name in SARAH: SSV): 

n(p 2 ) - cscccrY^F^p 2 , m 2 F ,m 2 ) (D.45) 

5. Vector boson loop (two 3-point interactions, generic name in SARAH: SVV): 

H(p 2 ) = csccCR-T^Bvip 2 ^^) (D.46) 

6. Vector boson loop (4-point interaction, generic name in SARAH: SSVV): 

Tl(p 2 ) ^ c s c c TA (m 2 v ) (D.47) 

Corrections to vector boson. 

1. Fermion loop (generic name in SARAH: FFV): 

n T (p 2 ) - c s c c c R {{\Y\\ 2 + \T R \ 2 )H (p 2 , m 2 v ,m 2 F ) + A^(TiT 2 R )B (p 2 , m 2 v , m 2 F )) (D.48) 

2. Scalar loop (generic name in SARAH: SSV): 

Il T (p 2 ) = -Acsc c c R \T\ 2 B 22 {p 2 ,m 2 Si ,m 2 2 ) (D.49) 

3. Vector boson loop (generic name in SARAH: VVV): 

U T (p 2 ) = \r\ 2 c s ccCR(-(4p 2 + m 2 Vl +m 2 v jB (p 2 ,m 2 Vl ,m 2 v J-SB 22 (p 2 ,m 2 Sl ,m 2 2 )) (D.50) 

4. Vector-scalar loop (generic name in SARAH: SVV): 

Il T (p 2 ) = \T\ 2 c s ccCRBo(p 2 ,m 2 v ,m 2 ) (D.51) 

We need here only the diagrams involving three-point interactions because the four-point interactions are 
related to them due to gauge invariance. 
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Appendix D.2.3. One-loop corrections to masses 

The one-loop self-energies can be used to calculate the one-loop masses and mass matrices. For the 
one-loop corrections of scalars, the radiatively corrected mass matrix is 



m lL (Pi) = m T b ~ U Ss(pf), 

while the one-loop mass of a vector boson V is given by 

m i'L(Q) = m T V + Re{lTy y (m^ v ')}. 



2,S 



(D.52) 



(D.53) 



According to the conventions of the counter terms of [l(|, the one-loop mass matrices M^ L of Majorana 
lected tc 



fermions are connected to the one-loop self-energies and tree-level mass matrix by 

1 



In the case of Dirac fermions, the one-loop corrected mass matrix is 

Mf L {p\) = Mf - E+(p») - E+(^)M| + - M| + S+( Pl 2 ) 

Appendix E. The minimal supersymmetric standard model 

Appendix E.l. Vector super fields 



SF 


Spin -i 


Spin 1 


SU(N) 


Coupling 


Name 


B 




B 


U(l) 


5i 


hypercharge 


W 


■V 


W 


SU(2) 


92 


left 


9 




9 


SU(3) 


93 


color 



Appendix E.2. Chiral superfields 



SF 


Spin 


Spin i 


q 


Q 


q 


i 


I 


I 


Hd 


H d 


Hd 




H u 


H u 


d 


d* 
a R 


d* 
a R 


u 


&h 


U *R 


e 


e R 


e R 



Generations 



(17(1) «) SU(2)® SU(3)) 



(£,2,3) 

(-i2,i) 

(-1,2,1) 
(1,2,1) 
(|,1,3) 

(-1,1,3) 
(1.1,1) 



Appendix E.3. Superpotential and Lagrangian 
Superpotential. 



(D.54) 



(D.55) 



W=Y u uqH u ~Y d dqH d - Y e elH d +fiH u H d 



(E.l) 
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Soft-breaking terms. 

-L SB .w = - H^H^B^ + HJH+B^ + H^d* R ia 5 a pdL,jpT d ^j - H d d* R ia 5 a pu L ,jpT d ^ 

+ H d^*R,i^L,jT e ^j - e Ri VL,jT e ,ij - H+u Ria 5 a pd L jpT u sj + H®u Ria 5 a pu L jpT Uti j 

-L SB ,4> = + ™ 2 Hd \H°\ 2 + m 2 H jH-\ 2 + m 2 H JH°\ 2 + m 2 H jH+\ 2 + d* L ^8 a pm 2 qtij d LM 
+ d R ^J a f)m 2 d , l0 dR,3i3 + e* Lij m 2 i:j e Li i + e R ^m 2 elj e R ^ + u Ljl3 8 al 3m 2 qtJ u L , la 

-Lsb.x = \ (a|Mx + M 2 \ 2 y i + MzXlt + h.c.) 

Appendix E.3.1. Gauge fixing terms 
Gauge fixing terms for gauge eigenstates . 



L 



GF 



2£g 



2£w ' 



duW 1 



Gauge fixing terms for mass eigenstates after EWSB . 

Lgf = - ^-^7 - T^-d^Sa - 27~( - A\m z iz + - ^-(iH^m w -^ w + d^W'^j 



^p""' 2£ G ^ 2i z 
Appendix E.4. Vacuum expectation values 



1 



= ^ d + iad + Vd> 



H °^v^ ^ u + i<Ju + Vu ^ 



Appendix E.5. Rotations of vector bosons and gauginos after EWSB 



W7 = — W~ 
lp ^2 p 



— W + 

V2 P 



w. 



2p 



-i^=W~ +i-^W+ 



W 3p = cos QwZp + sin Qwlp , B p = cos Qwlp — sin Q>wZ p 

1 ~ 1 ~ 



Aw 3 = W a 



Appendix E.6. Rotations in matter sector to mass eigenstates after EWSB 

In the following, Greek letters eti,(3i refer to color indices and Oi,pi to generations indices. 



1. Mass matrix for neutralinos, basis: ^A^, W°, H^, H^j 



/Mi 
M 2 

V \gwu -\giv u 



-\g\Vd 


\g\v u 


\givd 


-\giv u 













This matrix is diagonalized by N: 



Nm^N^ = mf a 



with 



*2 



t 2 



*2 



t 2 
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2. Mass matrix for charginos, basis: (w ,H d ^j/ (\V + ,H^ 

( M 2 -^92V U 



This matrix is diagonalized by U and V 



-^92V d II 



with 

tl t2 t2 

3. Mass matrix for leptons, basis: (ez, i0l ) / (e* Rpi ) 

m e = ( -^v d Y e>Pl0l ^ 
This matrix is diagonalized by {/£ and £7j^ 

U2*m e Ug = rof° 

with 

'/ .- = £ ^-i^J . - £ 

*2 *2 

4. Mass matrix for down-quarks, basis: (dx, i0iai ) / ^d^ iPlJ g 1 ^ 

™d = ( 7f«d^ai^i5 / d,pioi 

This matrix is diagonalized by Uf and E/jj 

U d L *m d U d / = mr 

with 

C^L.ia = £ U L*ji D L,ja , d R . ia = £ U Rlj D* Rja 



«2 «2 



5. Mass matrix for up-quarks, basis: (ti£, 0iai ) / (w^p^) 

This matrix is diagonalized by and U R 

U u L *m u U^ = m d : a 

with 

UL.ia = £ U^'jiULja , Wfl iia = £ U Rij U Rj 



t 2 
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6. Mass matrix for down-squarks, basis: {d Lt0l0ll /d Rt020l2 ^ , (d* LpilSi ,d Rp2 

1 3 

mil = ^S aif3l {l2{2 m l 0ipi +v 2 d J2Yla p Y d , a0l ) - (zgl + gf) ( - v\ + t%)s oiPl ) (E.26) 

a—l 

mi2 = -j=8 ail 3 2 (v d T dtP20l - v u ^*Y d , p20 ^j (E.27) 
1 3 

W22 = j2<W 2 ( 6 ( 2m lp20 2 + W d ^ i d* 02 a i d :P2 a) + ffl ( ~ + ^)<5 02 p 2 ) (E.28) 

a=l 

This matrix is diagonalized by Z 15 : 

Z D m\Z D * = mf| (E.29) 

with 

dh,ia = X Zfi*dja , d,R,ia = ^ %fi* dja (E.30) 

*2 *2 

7. Mass matrix for sneutrinos, basis: (fx, 0l ) / (^l^) 

™2=( sKp. + ^ + ^-^ + ^W ) (E.31) 
This matrix is diagonalized by Z v : 

Z v m\Z v * = mfl (E.32) 

with 

^, = E^ (E.33) 

t 2 

8. Mass matrix for up-squarks, basis: {u L . 0iai , u R , 02a2 ) / {u* L . piPl ,u* R . P2fj ^j 

1 3 
mn = ^S aiPl (l2(2ro* j0ipi + ^ Xr„* api y„, aol ) - ( - 3 ff i + ( - «S + ^)<W) (E.34) 

a=l 

mi2 = -^=S ai p 2 ( - v d ^*Y UtP20l + v u T UtP2Q ^j (E.35) 
1 3 

m 2 2 = Q$a 2 f3 2 (3vlJ2 Y u,o 2 a Y u, P 2a + 6rnl, P2 o 2 +5i(-«2 + «3)*o 2P2 ) (E.36) 
a=l 

This matrix is diagonalized by Z 17 : 

^m?Z^ = m*? (E.37) 

with 

WL, M = X Z ji'*Uja , U R . ia = ^ Z^*Uj a (E.38) 

t 2 t 2 

9. Mass matrix for sleptons, basis: {e L , 0l , eR, Q2 ) / (e* L pi , e* Rp2 ) 

1 3 

mn = g (4«3 X r e % Pl i;,a 0l + 8m? i0lPl + ( - , 92 2 + gf) (-v 2 u + ^)<J 0lPl ) (E.39) 

a— 1 

m 12 = -^=(v d T eiP20l - v u /j,*Y etP20l ^ (E.40) 
1 3 

m 22 = - (2^ ^ Y e ; 2 /e, P2 a + 4m^ ?202 + .9? ( - ^ + ^) (5„ 2P2 ) (E.41) 

a=l 
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This matrix is diagonalized by Z E : 

Z E m\Z E * = m*f 

with 



7 E,*~ 

*2 t 2 



10. Mass matrix for scalar Higgs, basis: {(j>d, 4>u) 

- \ (g\ + 5 2 2 ) «d«„ ra\ u + | M | 2 - | (g\ + .g 2 ) ( - 3v 2 u + vj 



™l = 



This matrix is diagonalized by Z H : 

Z» m \Z»*=m^ h 

with 

^ = z Z z " h ^ <t>u = Y, z % h J 

*2 t 2 

The mixing matrix can be parametrized by 

_ ( — sin a cos a 



cos a sin a 

11. Mass matrix for pseudo scalar Higgs, basis: (ad,a u ) 



2 



+ H 2 + i(.9 2 + .9 2 2 )(-^ + ^) 



^ { + | M | 2 - I (.g 2 + .g 2 ) ( - v\ + « 

This matrix is diagonalized by Z A : 

Z A m\ Z A * = mii a A0 

with 

The mixing matrix can be parametrized by 



m H- = 



sin /3 cos /3 
12. Mass matrix for charged Higgs, basis: 

+ M 2 + g + ffl) («3 - « 2 ) i3 2 2 ^«„ + B M 

\g\v d v u + B* m 2 ^ + | M | 2 + | (.g 2 - .g 2 ) (i; 2 - «; 

This matrix is diagonalized by Z + : 

Z+m 2 H -Z+J = m$%- 

with 

*2 *2 

The mixing matrix can be parametrized by 



^ + _ / — cos/3 sin/3 
y sin /3 cos /3 
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Appendix E. 7. Tadpole equations 



= | ( - 8 Wd K(X) + S^l^l 2 + ^(Sm^ - .g 2 ^ + g\vl - g 2 2 v 2 d + .gfw 2 )) (E.57) 

Appendix F. Particles and parameters of the MSSM in SARAH 

Particles 

We show here only the eigenstates after EWSB 
1. Fermions 





Cha [{generat ion}] 




Chi [{generation}] 


diet 


Fd [{generation, color}] 




Fe [{generation}] 




Fu [{generation, color}] 




Fv [{generation}] 


9i 


Glu [{generation}] 







2. Scalars 



3. Vector bosons 



di a Sd [{generation, color}] 
Ui a Su [{generation, color}] 
hi hh [{generation}] 
H~ Hpm [{generation}] 


i)i Sv [{generation}] 
e~i Se [{generation}] 
A" Ah [{generation}] 




gip VG [{generation, lorentz}] 
-y p VP [{lorentz}] 


W~ VWm [{lorentz}] 
Z p VZ [{lorentz}] 



4. Ghosts 



Vi 


gG [{generation}] 




gWm 


v + 


gWmC 


rfi 






gz 







Parameters 



9i 


gl 


92 


g2 


93 


g3 


Y u 


Yu 


T u 


T[Yu] 


Yd 


Yd 


T d 


T [Yd] 


Y e 


Ye 


T e 


T [Ye] 




\[Mu] 


B, 


B [\ [Mu] ] 


m\ 


mq2 


mf 


ml2 


m \ 


mHd2 


m H u 


mHu2 


m l 


md2 


m 2 u 


mu2 


2 

m e 


me2 


Mi 


MassB 


M 2 


MassWB 


M 3 


MassG 


v d 


vd 


v u 


vu 


Ow 


ThetaW 


4>g 


PhaseGlu 


Z D 


ZD 


Z v 


ZV 


Z U 


ZU 


Z E 


ZE 


Z H 


ZH 


Z A 


ZA 


z+ 


ZP 


N 


ZN 


u 


UM 


V 


UP 


u c L 


ZEL 




ZER 




ZDL 


u d R 


ZDR 




ZUL 




ZUR 


a 


\ [Alpha] 


p 


\ [Beta] 
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Appendix G. Changes in comparison to version 1 of SARAH 

We want shortly give an overview to the user about the main changes in the new version of SARAH. 

1. New physical output: 

(a) One- and two-loop renormalization group equations. 

(b) One-loop self-energies and one-loop corrected tadpoles. 

(c) All irreducible representations of chiral superfields possible . 

(d) Some representation theory with regard to SU(N) gauge groups. 

(e) Implicit charge indices are no longer restricted to SU(2)l. 

(f) Check for charge conservation of a model. 

2. Changes in definition of models: 

(a) DEFINITION statements to structure the model file. 

(b) Definition of global properties of parameters and particles. 

(c) Additional interactions added to the Lagrangian also rotated to new basis. 

(d) Possibility to add phases to particles. 

(e) Possibility to decompose one field with several generations in flavor eigenstates. 

(f) Soft-breaking terms are named to SLHA 2 conventions. 

(g) Improved routines to read LesHouches files. 

3. Changes in output 

(a) Speed of CalcHep/CompHep and I^TgXoutput significantly improved. 

(b) Adding of running coupling in CalcHep model file. 

(c) Possible suppression of splitting of four-scalar interactions in CalcHep. 

(d) Typesetting in IM^X output improved. 

(e) Overview of superfields, important terms of the Lagrangian, mass matrices, tadpole equations, 
RGEs and self-energies also added to the I^TEXoutput. 

(f) New commands MassMatrix and TadpoleEquation to have easier access to information. 

4. New models 

(a) MSSM with trilinear R-parity violation. 

(b) Singlet Extended Minimal Supersymmetric Standard Model (MSSM). 

(c) The [/(l)-Extended Minimal Supersymmetric Standard Model (UMSSM). 

(d) The Secluded [/(l)-Extended Minimal Supersymmetric Standard Model (sMSSM). 

(e) The near-to-Minimal Supersymmetric Standard Model (nMSSM). 
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